Continuing my expolartion of Notification of Infectious Disease data

https://journocode.com/2016/01/28/your-first-choropleth-map/

library("rgdal")
LocalAuthorities<-readOGR("https://opendata.arcgis.com/datasets/686603e943f948acaa13fb5d2b0f1275_4.geojson")
OGR data source with driver: GeoJSON 
Source: "https://opendata.arcgis.com/datasets/686603e943f948acaa13fb5d2b0f1275_4.geojson", layer: "686603e943f948acaa13fb5d2b0f1275_4"
with 380 features
It has 10 fields

“Statutory notifiable diseases - number of cases reported in week 51 of 2017”

library("reshape2")
library("magrittr")
geo_disease_data<-read.csv("geo_disease_data.csv")
geo_disease_data[is.na(geo_disease_data)] <- 0
geo_disease_data$LocalAuthority<-gsub(" UA", "", geo_disease_data$LocalAuthority)
LA_Diseases<-sp::merge(LocalAuthorities, geo_disease_data, by.x="lad16nm", by.y="LocalAuthority")
library("leaflet")
diseases <- names(geo_disease_data)[-1]
max_cases <- max(melt(geo_disease_data)$value)
Using LocalAuthority as id variables
pal <- colorNumeric(
   palette = "viridis",
   domain = c(0,max_cases)
 )
myleaflet<-leaflet(width = "100%") %>% 
  addProviderTiles("Esri.WorldGrayCanvas") %>%
  addLegend( pal = pal, 
             values = 0:max_cases,
             title = "Cases",
             opacity = 1) %>%
  addLayersControl(baseGroups=diseases, 
                   position = "bottomleft", 
                   options = layersControlOptions(collapsed = FALSE))
for (active_disease in diseases) 
  {
  myleaflet <- myleaflet %>%
    addPolygons(data=LA_Diseases, 
              fillColor=~pal(LA_Diseases[[active_disease]]),
              fillOpacity = 0.8,
              color = "black",
              weight = 1,
              popup = paste(LA_Diseases$lad16nm, LA_Diseases[[active_disease]],"cases"),
              group = active_disease
              )
  }    
myleaflet
Unique_LA_Map<-data.frame(LocalAuthority=unique(LocalAuthorities$lad16nm) %>% sort())
Unique_LA_Diseases<-data.frame(LocalAuthority=unique(geo_disease_data$LocalAuthority) %>% sort())
In_Map_Not_Diseases<-data.frame(LocalAuthority=
  Unique_LA_Map$LocalAuthority[!(Unique_LA_Map$LocalAuthority %in% Unique_LA_Diseases$LocalAuthority)])
In_Diseases_Not_Map<-data.frame(LocalAuthority=
  Unique_LA_Diseases$LocalAuthority[!(Unique_LA_Diseases$LocalAuthority %in% Unique_LA_Map$LocalAuthority)])
LS0tDQp0aXRsZTogIlIgTm90ZWJvb2siDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KDQpDb250aW51aW5nIG15IGV4cG9sYXJ0aW9uIG9mIE5vdGlmaWNhdGlvbiBvZiBJbmZlY3Rpb3VzIERpc2Vhc2UgZGF0YQ0KDQpodHRwczovL2pvdXJub2NvZGUuY29tLzIwMTYvMDEvMjgveW91ci1maXJzdC1jaG9yb3BsZXRoLW1hcC8NCg0KYGBge3J9DQoNCmxpYnJhcnkoInJnZGFsIikNCg0KTG9jYWxBdXRob3JpdGllczwtcmVhZE9HUigiaHR0cHM6Ly9vcGVuZGF0YS5hcmNnaXMuY29tL2RhdGFzZXRzLzY4NjYwM2U5NDNmOTQ4YWNhYTEzZmI1ZDJiMGYxMjc1XzQuZ2VvanNvbiIpDQoNCmBgYA0KIlN0YXR1dG9yeSBub3RpZmlhYmxlIGRpc2Vhc2VzICAtIG51bWJlciBvZiBjYXNlcyByZXBvcnRlZCBpbiB3ZWVrIDUxIG9mIDIwMTciCQkJCQkJCQkJDQoNCmBgYHtyfQ0KDQpsaWJyYXJ5KCJyZXNoYXBlMiIpDQpsaWJyYXJ5KCJtYWdyaXR0ciIpDQoNCmdlb19kaXNlYXNlX2RhdGE8LXJlYWQuY3N2KCJnZW9fZGlzZWFzZV9kYXRhLmNzdiIpDQoNCmdlb19kaXNlYXNlX2RhdGFbaXMubmEoZ2VvX2Rpc2Vhc2VfZGF0YSldIDwtIDANCg0KZ2VvX2Rpc2Vhc2VfZGF0YSRMb2NhbEF1dGhvcml0eTwtZ3N1YigiIFVBIiwgIiIsIGdlb19kaXNlYXNlX2RhdGEkTG9jYWxBdXRob3JpdHkpDQoNCkxBX0Rpc2Vhc2VzPC1zcDo6bWVyZ2UoTG9jYWxBdXRob3JpdGllcywgZ2VvX2Rpc2Vhc2VfZGF0YSwgYnkueD0ibGFkMTZubSIsIGJ5Lnk9IkxvY2FsQXV0aG9yaXR5IikNCg0KYGBgDQoNCmBgYHtyfQ0KDQpsaWJyYXJ5KCJsZWFmbGV0IikNCg0KZGlzZWFzZXMgPC0gbmFtZXMoZ2VvX2Rpc2Vhc2VfZGF0YSlbLTFdDQoNCm1heF9jYXNlcyA8LSBtYXgobWVsdChnZW9fZGlzZWFzZV9kYXRhKSR2YWx1ZSkNCg0KcGFsIDwtIGNvbG9yTnVtZXJpYygNCiAgIHBhbGV0dGUgPSAidmlyaWRpcyIsDQogICBkb21haW4gPSBjKDAsbWF4X2Nhc2VzKQ0KICkNCg0KbXlsZWFmbGV0PC1sZWFmbGV0KHdpZHRoID0gIjEwMCUiKSAlPiUgDQogIGFkZFByb3ZpZGVyVGlsZXMoIkVzcmkuV29ybGRHcmF5Q2FudmFzIikgJT4lDQogIGFkZExlZ2VuZCggcGFsID0gcGFsLCANCiAgICAgICAgICAgICB2YWx1ZXMgPSAwOm1heF9jYXNlcywNCiAgICAgICAgICAgICB0aXRsZSA9ICJDYXNlcyIsDQogICAgICAgICAgICAgb3BhY2l0eSA9IDEpICU+JQ0KICBhZGRMYXllcnNDb250cm9sKGJhc2VHcm91cHM9ZGlzZWFzZXMsIA0KICAgICAgICAgICAgICAgICAgIHBvc2l0aW9uID0gImJvdHRvbWxlZnQiLCANCiAgICAgICAgICAgICAgICAgICBvcHRpb25zID0gbGF5ZXJzQ29udHJvbE9wdGlvbnMoY29sbGFwc2VkID0gRkFMU0UpKQ0KDQpmb3IgKGFjdGl2ZV9kaXNlYXNlIGluIGRpc2Vhc2VzKSANCiAgew0KICBteWxlYWZsZXQgPC0gbXlsZWFmbGV0ICU+JQ0KICAgIGFkZFBvbHlnb25zKGRhdGE9TEFfRGlzZWFzZXMsIA0KICAgICAgICAgICAgICBmaWxsQ29sb3I9fnBhbChMQV9EaXNlYXNlc1tbYWN0aXZlX2Rpc2Vhc2VdXSksDQogICAgICAgICAgICAgIGZpbGxPcGFjaXR5ID0gMC44LA0KICAgICAgICAgICAgICBjb2xvciA9ICJibGFjayIsDQogICAgICAgICAgICAgIHdlaWdodCA9IDEsDQogICAgICAgICAgICAgIHBvcHVwID0gcGFzdGUoTEFfRGlzZWFzZXMkbGFkMTZubSwgTEFfRGlzZWFzZXNbW2FjdGl2ZV9kaXNlYXNlXV0sImNhc2VzIiksDQogICAgICAgICAgICAgIGdyb3VwID0gYWN0aXZlX2Rpc2Vhc2UNCiAgICAgICAgICAgICAgKQ0KICB9ICAgIA0KDQoNCm15bGVhZmxldA0KDQpgYGANCg0KYGBge3J9DQpVbmlxdWVfTEFfTWFwPC1kYXRhLmZyYW1lKExvY2FsQXV0aG9yaXR5PXVuaXF1ZShMb2NhbEF1dGhvcml0aWVzJGxhZDE2bm0pICU+JSBzb3J0KCkpDQpVbmlxdWVfTEFfRGlzZWFzZXM8LWRhdGEuZnJhbWUoTG9jYWxBdXRob3JpdHk9dW5pcXVlKGdlb19kaXNlYXNlX2RhdGEkTG9jYWxBdXRob3JpdHkpICU+JSBzb3J0KCkpDQoNCkluX01hcF9Ob3RfRGlzZWFzZXM8LWRhdGEuZnJhbWUoTG9jYWxBdXRob3JpdHk9DQogIFVuaXF1ZV9MQV9NYXAkTG9jYWxBdXRob3JpdHlbIShVbmlxdWVfTEFfTWFwJExvY2FsQXV0aG9yaXR5ICVpbiUgVW5pcXVlX0xBX0Rpc2Vhc2VzJExvY2FsQXV0aG9yaXR5KV0pDQoNCkluX0Rpc2Vhc2VzX05vdF9NYXA8LWRhdGEuZnJhbWUoTG9jYWxBdXRob3JpdHk9DQogIFVuaXF1ZV9MQV9EaXNlYXNlcyRMb2NhbEF1dGhvcml0eVshKFVuaXF1ZV9MQV9EaXNlYXNlcyRMb2NhbEF1dGhvcml0eSAlaW4lIFVuaXF1ZV9MQV9NYXAkTG9jYWxBdXRob3JpdHkpXSkNCg0KYGBgDQoNCg==